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Experimental tests of Bell-type inequalities distinguishing between quantnm mechanics and local 
realistic theories remain of considerable interest if performed on massive particles, for which no 
conclusive result has yet been obtained. Only two-particle experiments may specifically test the 
concept of spatial nonlocality in quantum theory, whereas single-particle experiments may generally 
test the concept of quantum noncontextuality. Here we have performed the first Bell-type experiment 
with a beam of thermal-neutron pairs in the singlet state of spin, as originally suggested by J.S. Bell. 
These measurements confirm the quantum-theoretical predictions, in agreement with the results of 
the well-known polarization experiments carried out on optical photons years ago. 
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Since the starting point of Bell’s theorem on non¬ 
locality of the quantum theory in 1965 [T], theoretical 
speculations and search for experimental verifications of 
various inequalities derived from the theorem have been 
pursued over the years. Most of the performed experi¬ 
ments were based on correlations between the polariza¬ 
tion of pairs of photons, in particular high-energy pho¬ 
tons produced by positronium annihilation [5] , or optical 
photons emitted by atomic cascade [31 H or by para¬ 
metric down-conversion [S]. One test was also done on 
massive particles to measure the spin correlation in pro¬ 
ton pairs prepared in low energy proton-proton scattering 
[5] . These early experiments on positronium annihilation 
and proton-proton scattering did not give conclusive re¬ 
sults, mainly because of the lack of efficient linear polar¬ 
ization filters both for protons and high-energy photons, 
and because of other difficulties related to the apparatus. 
One more recent and complex test performed on massive 
particles, pairs of entangled ^Be"*" ions, appears as con¬ 
ceptually questionable because the manipulations of the 
two particles of a pair were not done individually and 
independently of each other, as would be required to ver¬ 
ify spatial non locality [7]. The extensive tests working 
with atomic cascades producing pairs of optical photons 
have been the most successful (together with those com¬ 
pleted afterward by using parametric down-conversion), 
mainly because of the availability of very efficient polar¬ 
ization filters and an efhcient technique for the manipula¬ 
tion of visible light beams. The experimental limitations 
of such experiments are mainly related to the relatively 
low quantum efficiency of the detectors and the difficulty 
to select and detect both members of a photon pair, be¬ 
cause the variable recoil momentum of the emitting atom 
diminishes their direction correlation. However, with a 
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supplementary assumption justifying that the detected 
photon pairs constitute a statistically representative sam¬ 
ple of the whole ensemble of photons emitted from the 
source, these experiments confirm the predictions of the 
quantum theory with great accuracy. A detailed and 
comprehensive analysis of the three classes of test ex¬ 
periments performed on Bell’s theorem (e.g. atomic cas¬ 
cades, proton-proton scattering, positronium decay), is 
given in [5]. 

The above considerations, and particularly the limits 
of the experiments with massive particles, have moti¬ 
vated us to perform a Bell-type experiment working with 
a collimated beam of thermal neutrons prepared in such a 
way that, within the experimental limitations of the mea¬ 
surement, neutron pairs in the singlet state of spin were 
identified by the detecting system. We recall that an ex¬ 
perimental test on such a physical system is particularly 
interesting since Bell’s inequality is violated by a large 
amount in the singlet state. We also recall that beau¬ 
tiful interferometric measurements have been done with 
single-neutron beams, testing correlations between two 
degrees of freedom of the single particle [MI]. However 
such experiments are not equivalent to those performed 
with two-neutron beams, which refer to the issue of spa¬ 
tial nonlocality for two separated particles. This latter 
concept is meaningless for a single particle, and the in¬ 
terferometric experiments have rather tested the general 
concept of quantum noncontextuality. 

Following the results of the previous measurements 
on the fermion anti-bunching [m Ha], the present 
measurements have been carried out at the Institute 
Laue Langevin (Grenoble, France) by using the pri¬ 
mary spectrometer of instrument INIO, which produces 
a monochromatic beam of thermal neutrons from an al¬ 
most perfect Si(lll) single crystal in the almost perfect 
backscattering configuration. The neutron intensity at 
the test position was 2 x 10^ n/s on a beam size of 3 x 4 
cm^. The energy spread cannot be measured directly but 
it can be estimated from the monochromator geometry 



2 


to be AE < 0.02 ^eV, corresponding to a coherence time 
Tc = h/AE > 20 ns at 6.27 A. The distance of the mul¬ 
tianode detector from the monochromator was ^ 12 m 
and ^ 2 m from the collimator exit; the whole equipment 
was within the shielding of INIO. 

The extraction of a beam of correlated particles in the 
singlet state can be described as follows. The nature of 
the emission of thermal neutrons in the source (the neu¬ 
tron moderator of ILL) is Poissonian, so that there is a 
small but finite probability of having two neutrons within 
the detection time of our apparatus. The beam emitted 
from the source is monochromatized by reflection on a 
quasi-perfect Si crystal. The beam being spin unpolar¬ 
ized, it is equally likely that a neutron pair emitted within 
the coherence time t,. of the monochromator will either 
occur in one of the three triplet states or in the singlet 
state, i.e. the triplet states will occur 3/4 of the time. 
In a gas of fermions there is the tendency for particles 
of the same spin to avoid each other, a tendency arising 
from the exchange antisymmetry of the wave function. 
More specifically the two neutrons of a pair emitted in 
the triplet state and traveling along the long collimators 
(^10 m) from the monochromator to the polarizers are 
in the same spatial mode (antisymmetrical spatial wave- 
function), so that they avoid each other scattering in di¬ 
rections different from the one of the emerging collimated 
beam, and they will not reach the detectors. The neu¬ 
trons of a pair in the singlet state do not avoid each other 
traveling along the collimators; they will pass through the 
two separated polarizers, and after them, are distinguish¬ 
able particles that will be detected without exhibiting 
interference effects. Consequently, within a time inter¬ 
val of the order of Tc, only neutron pairs emitted from 
the source in the singlet state can be detected in our ap¬ 
paratus. This is the physical effect already measured in 
[HIS], and applied in the present experiment. 

The conceptual scheme of the experiment was similar 
to the one described in | 1 ^ . with the addition of two 
polarizers (1 and 2 in FigQ, each one mounted along 
two slightly different neutron paths from the beam out¬ 
put window. Both polarizers consisted of a disk, 3.8 cm 
diameter and 2.6 cm thick, made of FesAl intermetallic 
compound and mounted between the poles of a 0.2 T per¬ 
manent magnet. The FeaAl inter-metallic compound was 
chosen because the first Bragg edge, corresponding to the 
( 111 ) reflection of its DO 3 structure, is at A = 6.67 A, 
while the second Bragg edge due to the (200) reflection 
is at A = 5.78 A. So, at the working neutron wavelength 
of 6.67 A the transmission of the disk is due only to the 
capture cross section of Fe and to the (111) Bragg scat¬ 
tering. If a polarizer is magnetically saturated, the cross 
section is large when the Sz component of the neutron 
spin is parallel to the magnetization, and it is small when 
Sz is antiparallel M- With our 0.2 T magnet the mea¬ 
sured transmission coefficients were = 22 . 1 % and e”l= 
8.4%. In the present apparatus, polarizer 2 was mounted 
on a motorized rotation stage and could be rotated with 
respect to polarizer 1 mounted within a fixed horizontal 



Figure 1: Scheme of the experiment. Distances along the 2 - 
direction are not on scale. Sizes of the neutron guide and 
collimating window are respectively 3x3 cm^ and 3x4 cm^. 
64CH = 64 channels from the multi-anode pixels; DAQ = 
data acquisition system; OR-module is a logical circuit trig¬ 
gering DAQ on receiving a signal from whatever channel of 
the 64CH discriminator module. The polarizers are separated 
vertically, fronting the two halves of the detectors; erroneous 
visual impression may arise from the drawing perspective. 


field. The relative angle 'd of their polarization axes was 
variable in the range from 0° to 360° ± 0.1°. Since the 
polarizers can be considered as state preparation devices, 
their position along the neutron path was not important, 
provided that the two neutron paths were identical as it 
was in our apparatus. 

As in [13] , the position-sensitive detector was a Hama¬ 
matsu H8500 multi-anode photomultiplier having 8 x 8 
anode-pixels with a pixel size of 5.8x5 .8 mm^. The 
photomultiplier was coupled to a 0.15 mm thick lithium 
glass scintillator (®Li 98% enriched) directly coupled to 
the anode window with a thin layer of optical grease. 
The scintillator was divided into two different parts with 
a separation distance of ~ 0.4 mm, so that each part 
of the scintillator was looking at the neutrons passing 
through either polarizer 1 or polarizer 2 , and there was 
no cross-talk between the two detecting areas due to the 
neutrons impinging on the separation area. Specific pixel 
groupings could be defined during the off-line analysis to 
optimize the view of the upper and lower parts of the 
neutron beam, so that each of the two detectors, Di and 
D 2 j recorded the neutrons passing through one specific 
polarizer. Generally each detector was formed of a row 
of 4-6 pixels with total area < 3.0x 0.5 cm^, chosen so 
as they were illuminated only by those neutrons that, fil¬ 
tered respectively by polarizer 1 and polarizer 2 , could 
not be intercepted by the Fe magnetic poles of polar¬ 
izer 2 at any angle d. The neutron arrival time at each 
pixel was measured by using an internal 40 MHz clock 
which provided 25 ns time resolution, shorter than both 
the light decay time of the lithium glass scintillator (~ 
250 ns) and the neutron total travelling time through 
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the scintillating sheet (~ 240 ns). The data acquisition 
system monitored events during repeated cycles of 10 s 
each, with an effective duty cycle of about 99.90%, where 
the dead time due to the VME read out cycle was of the 
order of 10 ms. 

In order to determine the spin correlation of two neu¬ 
trons, a time stamp was attributed to the neutron signals 
from the two detectors, so that the neutron pair was time 
correlated when two neutrons arrived at the two detectors 
Di and D 2 with time-stamp ti =start and t 2 =stop. In 
the off-line analysis we could study the coincidence rate 
as function of d = ^2 ~ ^i, or in other words as a function 
of the virtual spatial separation Uther<5 along the propaga¬ 
tion direction z. In [Ellla] it is also shown that, thanks 
to the specific characteristics of the INIO primary spec¬ 
trometers, there is a good neutron correlation when the 
two detector Di and D 2 are separated in the transverse 
direction up to a couple of centimeters. The issue of the 
existence or amplitude of macroscopic lateral coherence 
and its possible effects on coincidence measurements is 
controversial. For this reason we prefer to limit ourselves 
to report our measured value. We have no conclusive in¬ 
terpretation of this results, and we feel that this effect 
requires further experimental investigation. 

According to the above description, the experimental 
results that we are presenting provide the number of neu¬ 
trons monitored by £>2 at a time <2 delayed by 6 with 
respect to time ti when the first neutron is detected at 
Di. It is interesting to obtain a calculation of this number 
from the convolution of the predicted neutron-correlation 
function c(d,(5), and the response function of the detect¬ 
ing system w{5). 

The function w(5) = W exp[—describes the time 
broadening of the detecting system due to the scintilla¬ 
tor thickness and to the decay time of the emitted light. 
Actually, from the capture probability profile of the neu¬ 
tron along the scintillator thickness, the average capture 
time is about tc = 80 ns with an additional broadening 
from the light decay curve of the order of = 250 ns. By 
considering that both effects can be roughly described by 
an exponential function, the total decay time is expected 
to be tu, = [^/tc + = 60 ns. The choice of the 

Gaussian form of c{'d, h) = 1 — a('d) exp[—(5^/(2 t,?)] was 
dictated first by the need of describing reasonably (sim¬ 
ilarly to what was done in [T^]) the antibunching effect 
which produces, for 5 = 0 and for any value of d, a 3/4 
depression of the coincidences counts with respect to the 
random coincidences occurring for 6 ^ Tc] and second 
by the need to add also the simultaneous depression pro¬ 
duced by the presence of the polarizers on the neutrons 
pairs in the singlet state (see prediction Q, recalled later 
in the text). The value of the coefficient a{'d) must be 
3/4-I-I/4 = 1 for d = 0°. From the convolution [rt;*c]c(t') 
of w and c, the number of coincidences is then given by: 
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Figure 2: G(i9 = 0,(5) is shown as a function of the neutron 
time-delay S = t 2 —ti along the longitudinal z-direction; full 
square = experimental points. Two 4-pixel rows with relative 
separation 0.4 cm, A=150 ns. 



where C'('d, S) is the number of experimental coincidences 
as a function of 6, C(i?,oo) is the number of experimen¬ 
tal coincidences for h 3 > Tc, A is the coincidence time- 
window of the order of Tc, and W = 1. In Fig. the 
experimental results are compared with the prediction of 
Eq. 0. Tests taking a(0), and Tc as free parameters 
confirmed our estimates, and we found q:(0) = 1 ± 0.1, 
= 60 ± 9 ns and Tc = 78 ± 10 ns, in agreement with 
the previous estimates [Oj . 

We can now define the function: 


£12 (d) 


g(d,0) 

C{d,oo) 


( 2 ) 


which is independent of all possible fluctuations of the 
beam intensity and other sources of instrumental uncer¬ 
tainties. It represents the experimental value of the prob¬ 
ability of detecting the simultaneous occurrence, after 
transmission through the linear polarization analyzers, 
of the two neutrons of a pair emitted in the singlet state, 
with respect to all four states of the physical system. The 
experimental results obtained for Pi 2 (d) at five different 
values of d are shown in Fig. The acquisition time of 
each experimental point was in the range 18-24 h. 

In order to obtain the quantum-theoretical prediction 
for the outcome of the experiment, let us first calculate 
the probability £ 12 ( 1 ?) for perfect analyzers that transmit 
only neutrons with spin antiparallel to the magnetic field. 
During a given time interval, the source emits N neutron 
pairs. For N large enough, we define: Ni/N = Ci, W 
being the number of counts at detector-z; and we also 
recall that W = N/2 is the number of neutrons trans¬ 
mitted through perfect polarizer i when S = t 2 — ti ^ Tc- 
Since the fraction of the neutron pairs emitted in the sin¬ 
glet state and recorded after the polarizers at short time 
delays is 1/4 of all the pairs, we obtain 
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= = ( 3 ) 

where |'0s) is the singlet wave function and Qj is the 
operator which projects the spin of particle (i) on the di¬ 
rection of magnetic field with value -h/2. Yet a more 
realistic analysis, taking account of the limitations of 
the apparatus, can be adopted for comparison with the 
real experiment. We can adapt to result the analysis 
adopted by Clauser, Horne and Shimony (CHS) [SI HI] 
for pairs of optical photons with parallel linear polar¬ 
ization. They derive the following prediction: Pi 2 {'&) = 
ciC 2 [s^e 2 + 6^62 cos(2i?)]/4 and Pi 2 (oo,oo) = C 1 C 2 with 
-I- £™, e“ = — e™, where and e™ are 

respectively the maximum and minimum transmittance 
of polarizer i. The symbol 00 defines the absence of a 
polarizer. 

By adapting their approach, particularly the one in 
[8|, to our system of thermal neutrons, we have not per¬ 
formed a direct mathematical derivation of their result, 
because the significative differences to be considered be¬ 
tween the two cases have dictated to introduce changes 
in their calculation. First we recall that a linearly polar¬ 
ized photon propagating along the z axis has polarization 
direction in the x-y plane; differently a neutron propagat¬ 
ing along the z axis has spin components different from 
0 only in one x-y semiplane. We also recall that the op¬ 
tical linear polarizers have the maximum transmittance 
in direction perpendicular to that of minimum transmit¬ 
tance, while thermal neutron polarizers have maximum 
and minimum transmittance in antiparallel directions. 
Second, only coincidences occurring within a time inter¬ 
val S < Tc are considered in our analysis, i.e coincidences 
due to pairs of neutrons in the singlet state. 

We must now remark that a neutron polarizer like ours 
is a linear polarizer with oriented polarization axis, so 
that two of such perfect polarizers (ej = 1 and = 0) 
with parallel orientation of their axes cannot transmit 
two correlated neutrons in the singlet state (antiparal¬ 
lel spins), and coincidences may be observed only in the 
angular range 90° < i? < 180°. Consequently, the ideal 
result ^ can only be confronted with Eq. obtained 
later in the text, in this angular range. However, with 
real polarizers (sj = 0.221 and e| = 0.084 in our experi¬ 
ment) with oriented polarization axis, coincidence counts 
may also occur for 0° < < 90°, but still the probability 

of their occurrence cannot be confronted with prediction 

(§. 

As a consequence of such differencies between optical 
photons and neutrons, in the above CHS equation ef 
and e~ have become respectively ef and ej in the term 
proportional to cos('d). In addition, note that the neutron 
correlation function sin^(il/2)/2 [see Eq. §] substitutes 
the photon correlation function cos^{'d)/2; the symbol 00 
does not denote the absence of a polarizer, it denotes the 
absence of correlation between to two spins of a pair for 
S ^ Tc instead; the beam of neutrons in the singlet state 



Figure 3: Pi 2 ('d) as a function of the angle between the 
polarization axes. Two 4-pixel rows with relative separation 
0.4 cm, and A = 150 ns. Full line: quantum prediction (Eq. 

(§) 


is the fraction 1/4 of the whole ensemble emitted from 
the monochromator. 

Following the above considerations, we obtain for 0° < 
t? < 90° 


C(tl, 0 ) = 
C'('d, 00 ) 


t 


C 1 C 2 

4 


CiC2£‘i4 


and for 90° < < 180° 


(£^4 + 

4 


( 4 ) 


C(t?, 0 ) = 
C'('d, 00 ) 


C 1 C 2 

4 


S-t Sn 


ciC2s{el 

4 


(£41 -f£t£4cOs(t?) ' 

4 


( 5 ) 


where £*= £j-|-£j'. Following the definition of Eq. ([^ and 
assuming identical polarizers (£^2 = 4 ^ 2)1 derive 
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The results obtained by using Eq. Q are shown in Fig. 

(full line) as a function of'd. Clearly the theoretical 
prediction of the quantum theory is in agreement with 
the experimental points, within the errors of the present 
experiment. 

Recall also that for an analogous experiment which 
might be performed with ideal instrumentation and with 
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a source emitting only neutron pairs in the singlet state, 
the quantum prediction substituting Eq. ®> is : 

= = (7) 

and that this simple sinusoidal form of the ideal predic¬ 
tion (as other simple sinusoidal forms of quantum predic¬ 
tions, in other cases) is the sole origin for the violations 
of the CH inequality [T6| . 

Let us summarize the most significant experimental 
features of the present measurement: (1) from the highly 
monochromatic and collimated beam of thermal neutrons 
produced at the INIO beam of the ILL, by utilizing the 
anti-bunching effect of the neutrons in the triplet states, 
the two members of a pair in the singlet state could be 
detected at small time separation as a function of the 
angle between the two polarizers; (2) the effects on 
the data analysis due to the presence of undesirable and 
uncertain instrumental origin were greatly reduced by 
measuring Pi 2 ('i?) from the ratio of the coincidence rate 
at small time separation to the coincidence rate at large 
time separation {t 2 — ^ Tc). 

The measurement seems adequate to confirm the 
quantum-theoretical prediction ([^ which takes account 
of the real limitations of the present experiment, and, by 


extrapolation, the theoretical prediction Q for an ideal 
experiment performed with perfect instrumentation. 

We wish to remark explicitly that the efficiency of our 
polarizers was insufficient for a decisive confrontation 
of the experimental data with the Clauser-Horne (CH) 
inequality, confrontation which requires very high effi¬ 
ciency of all of the components of the apparatus. We 
also note that such limitation prevent our from perform¬ 
ing loophole-free measurements. On the other hand, we 
consider our result as a helpful contribution to this is¬ 
sue. Actually, the confirmation of quantum prediction 
(§ could but be reinforced by experiments performed 
with greater efficiency and statistical precision, as it was 
shown in the optical polarization experiments, and also 
commented by Bell Hi!: ”.. it is hard for me to believe 
that quantum mechanics works so nicely for inefficient 
practical set-up and is yet going to fail badly when suf¬ 
ficient refinements are made ...” Therefore, our results 
seem adequate to verify the predictions of the quantum 
theory for massive particles in the spin singlet state, and, 
although by extrapolation from the real measurement, to 
confirm violations of the CH inequality. 

We are grateful to P. Facchi and S. Pascazio for enlight¬ 
ening discussions and we acknowledge the Institut Laue 
Langevin of Grenoble for the beam time, which enabled 
the realization of the experiment. 
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